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We propose an experiment for the measurement of gravitational effect on cold atoms 
by applying a one-dimensional vertically sinusoidal oscillation to the magneto-optical 
trap; and observe the signature of low quantum energy shift of quantum bound states 
as a consequence of gravitational fluctuation. To this end, we present brief details of 
the experiment on a BEC, and a simplistic calculation of the Gross-Pitaevskii solution 
using Thomas- Fermi approximation with focus on the density of the BEC, for the time- 
dependent perturbation. Moreover, we calculate the power induced by quantum gravity 
on a generic atomic ensemble. We also address the possible challenges for the measure- 
ment of the expected results. And finally, we discuss the prospect of further developing 
this experiment towards measuring the effect of quantum spacetime fluctuations on cold 
atoms. 
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1. Introduction 

The use of cold atoms allows new ways of testing fundamental quantum phenomena 
( M.H. Anderson|1995 1. Among them is the elusive quantum gravity and its effect, which 



for years has avoided been understood and detected. Despite all the unsuccessful at- 
tempts, there have been suggestion that the low energy signature of quantum gravity 



may be detectable (C.H.-T. Wang 2012). According to the authors; and inspired by the 



success from the Lamb shift measurement in QED (see JR||1955 Bethe &: Brown|[l950 ), 



they believe there is a possibility of observing the energy shifting of the quantum bound 
state coupled to a quantum field under vacuum fluctuation; similar to that discovered 
by Lamb. Incidentally, due to the fact that the lowest energy of quantum gravity is the 
zero-point energy of vacuum, then by manipulating cold atoms or a BEC center of mass, 
formed in the lowest state, the energy shift may be amplified to the point of observation. 

In this paper we will present an experiment to investigate the behaviour of the BEC 
under such conditions, as required by the gravitational Lamb shift. Initially, a brief 
introduction about the gravitational Lamb shift mentioned above will be presented; then, 
we will describe the experiment to be cared out in the near future; and after, a brief 
analysis of the BEC behaviour and calculations of the Gross-Pitaevskii equation solution 
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under the harmonic perturbation focused on the density of the BEC. Moreover, we will 
calculate the general quantum mechanical case of an ensemble of atoms under the same 
harmonic perturbation. From there, we were able to observed that the ensemble atoms 
does indeed absorbs energy from the gravitational induced perturbation. Finally, we 
present the technical difficulties that may compromise the experiment in general, and 
also underline the results as well as the uncertainties of the experiment in the discussion. 



2. Gravitational Lamb Shift 

The idea of a gravitational Lamb shift according to (C.H.-T. Wang 2012) is based on 
the same principle used to determine the QED Lamb shift used for the Hydrogen atom. 
The paper used similar semiclassical derivation of the Lamb shift as derived by Wclton 



theory (M.O.Scully & Zubairy 19971. From such derivation they found that the energy 



shift of the condensate in a harmonic potential would be zero, unless the center of mass 
of the BEC is displaced from its original position. In that case, they found that the 
expected potential energy V is given by 



(AV) 
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(2.1) 



with m,mp,L being the atomic mass, Planck mass, and trap vertical length. Moreover, 
if they were using a frequency of lu = 2kHz with the number of atoms in the condensate 
equal N = 10 6 , then the energy ratio AE/hui would have increased to 0.5 %. Detailed 
information on their derivation of the gravitational Lamb shift is found in (see C.H.- 
T. Wang|[2012 ). 



3. Proposed experiment setup 

In this paper we propose an experiment for measuring quantum gravity based on cold 
atoms. To obtain a good guidance of the phenomenon, we will consider a simplified one- 
dimensional harmonic oscillation for the cold atoms and BEC. The trap used is a normal 
anisotropic spherical shape MOT. There is an illustration of the main equipments used in 
the experiment on fig.(l). The aim of the experiment is to achieve a reasonable displace- 
ment of the atom cloud or BEC center of mass from its origin, with a specific frequency. 
Such displacement, in our set up, provides a mechanical vertical sinusoidal oscillation of 
the MOT assemble due to an electro-acoustic vibration generator. 

The trap displacement is up to 1 cm (susceptible to change), giving the sine wave an 
amplitude of the same value. Incidentally, we assume that all the necessary conditions 
for the creation of BEC are met as well as the structure of the anisotropic harmonic trap 
is maintained under constant oscillations, that is, constant at the origin and during the 
oscillation of the system. The sinusoidal varying system will allow the formation of a 
harmonic perturbation of the metric due to the oscillating potential, which corresponds 
to a function of space and time in the direction of coordinate z given by 

z(t) = asin(fit) (3.1) 

where a is the amplitude and = 2irf the angular frequency of the oscillation. 
The acceleration is given by 



A(t) = -tt 2 asm(nt). 



(3.2) 
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Figure 1. Experiment setup, 1- trap, 2- quarter wave plate, 3- fibre optics, 4- tuneable laser, 5- 
electro-acoustic vibrator drive, 6-electro-acoustic vabration generator, 7- TOP maginetic coils, 
8- flexible vacuum pipe, 9- mirror, 10- ion pump, 11- trap magnetic coils (MOT), 12- telescope 



This gives rise to the time dependent potential perturbation as follows: 

V(z,t) = mA(t)z (3.3) 

using the BEC or atomic ensemble mass m. For this particular experiment we assume the 
phase to be zero and the frequency on the order of 1-2 kHz. The atom vapour species that 
we will be using is the isotope of Rubidium Rb87. The initial number of atoms in vapour 
will be of order 10 12 with the condensate having about 2000 to 10 6 atoms. In this case 
the interaction will be repulsive, that is, s-wave scattering length a s is positive. Within 
our 2.5 cm 3 trap, we would expect 2 x 10 3 - 2 x 10 4 atoms after the spatial distribution 
forming the isotropic thermal equilibrium condensate in the ground state (BEC). Then, 
by increasing the number of atoms in the trap to form a larger BEC, we will be able 
to observe if the corresponding effect of gravitational lamb shift will be proportional to 
that predicted by the theory. 

The 3D harmonic potential (i.e. TOP time orbiting potential of 7.5 kHz) generated by 
the magnetic coils are meant to be relatively symmetrically constant during the system 
oscillation. The time-dependent perturbation created by the vibration generator would 
allow the necessary conditions for the gravitational effect predicted to be exposed and 
detected. This vibration generator provides 5N of vertical force, and an oscillating fre- 
quency that varies from 1Hz to 15kHz, once controlled by the drive. 

Now, from the schematic above, it is worth mention that we did not show the RF 
bias magnetic field, external back side window of CCD camera, the front TOP magnetics 
coils , control systems, detection and the remaining Ioffee-Pritchard configuration mech- 
anisms. Incidentally, regarding the detection mechanism, we will be using a probing laser 
and a CCD camera to measure the level of energy, density and the decoherence of the 
condensate or atomic cloud. 
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4. Harmonic Perturbation on BEC 



We used the Gross-Pitaevskii equation (4.1 ) to study the behaviour of the BEC under 
the Harmonic perturbation, 

%n ^A = V 2 HXjt) + V ( X ,t)TiKx,t) + 4 ^W(x,t)\ 2 Tp(x,t). (4.i) 

There we have used V(X, t)x as the sum of the potential V(z, t) to the Anisotropic 
potential V(X, t) trap , and X = x,y, z coordinates: 



V(X, t) T = ^m a ojlx 2 + ]^n a uj 2 y y 2 + -m a uj 2 z z 2 + mA(t)z. (4.2) 



For a simplified analyse, we consider the numerical solution of the equation (4.1) using 
Thomas-Fermi Approximation. For this purpose, it is necessary that V(X, t)x is less than 
the chemical potential fi, in this paper we neglect the trap frequencies (i.e. the radial 
and axial frequencies uj x = ui y , and w z ). Here, m a is the atomic mass of Rb87, m the 
total mass of the BEC which will be Nom a , with Nq being the number of atoms in zero 
state, N the number of atoms. We also use an amplitude for the mechanical oscillation 
of 1cm for this particular case; the s-wave scattering length a s — 5.28nm, as for the Rb 
87. Incidentally, an observation of the effect of our potential using CPE solution towards 
the density of the BEC was obtained. Here, it was noticed that our potential generates 
a variational solution of GPE, and was acquired using Thomas-Fermi Approximation. 
Fig. (2), illustrated the behaviour of the matter wave density, as a function of time. 
There, it was noticed that as the mechanical oscillation evolved, the ratio of the density 
of Thomas- Fermi density tptf with n increased quite rapidly, and then collapsed to an 
imaginary solution of the wave function continuously. The Thomas- Fermi density Utf is 
given by( |Pethick fc Smith||2002[ ) 

n{z)=n TF = (fi-V[z,t))/g (4.3) 



where g = 4Trh 2 a s /m, fj, > V(X, t), with a s = 5nm, and mass of Rb87 BEC, /i ~ ng 
with n = /z(0)/fcsT c « 0.3 ( Dalfovo|1999 1 . The mechanical perturbation creates a rising 
effect on the total density of the condensate. We noticed that the effect forces the BEC to 
go beyond the Thomas-Fermi regime, giving negative density and therefore, forcing the 
solution of the wave function to be complex. From the graphic on the left of fig. (2), the 
density ratio increases to 0.5% above 1 with Ao = 1000, amplitude 1 cm and frequency 
1kHz; Whereas, for the graphic on the right, with the difference only in N — 10 6 , we 
observed and increase of 50%. These particular increases in the density ratio indicates 
that the Harmonic mechanical perturbation does affect the energy level to certain extent. 
However, a thorough analyse with different methods such as the Bogoliubov theory may 
be consider in future papers. 



5. Gravitational Power Absorption 

Here, we start by considering a one-dimensional harmonic oscillator in an initial state 
The harmonic oscillator is then perturbed by a small time-dependent potential 
V(x,t) = mA(t)x with a general (broadband) time dependent acceleration A{t). The 
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Figure 2. Density ratio variation with 1kHz, 1 cm amplitude, iVo = 1000 (left), and iVo = 10 6 

(Right) 



probability of finding the oscillator in the state |n) at a later time can be found using 
perturbation theory. 

We therefore consider the time-dependent Hamiltonian given as, 



H = H + V(x,t) 



(5.1) 



where 



V(x,t) = mA(t)x. 



(5.2) 



The generic time-dependent Hamiltonian in equation (5.1 ) does not have a typical ana- 



lytical solution. Using perturbative analysis, the basis coefficients c„(t) can be expanded 
in powers of the interaction. 



C n (t) = c^ + cW(t)+c^(t) + 



(5.3) 



For the generic potential, different cases of A(t) can be considered. For example, a 
harmonic potential can be given by 



A(t) = A sin(J2£). 



(5.4) 



It follows that 



V ni (t') = m{n\V(x,t)\i). 



(5.5) 



Substituting equation (5.2 1 into the above gives 



V n i(t') = mA(t) (n\x\i) . 



(5.6) 



Since 
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(n\x\i) 



2mw 

h 

2mw 



(n\a + a'\i) 

(y/i + l<J„,j+i + Vi5 n;i -i) 



we can calculate the first order perturbation coefficient as follows: 

c2 ) (t) = -^*dt / e^ t V ni (0 



Substituting equation ( |5.6[ ) into this equation yields 



■{n\x\i)A t {u n i) 



(5.7) 



(5. 



(5.9) 



where A t {uj n i) is the Fourier transform of the complex stationary stochastic function 
A(t), given by 



A t (u ni ) = / dt'e^'A{t'). 



The mean power spectral density of this integral is normally given by 



1 



(At(uj n i)A t (uJni)) 



(5.10) 



(5.11) 



where (A*(w n i)A t (ui n i)) denotes the statistical average. Additionally, It is known that 

(5.12) 



S(u) n i) = lim S t (uJ n i) 
t— too 

And the transition probability can be found using (|5.9|) to be 



Pni — IC, 



(1)|2 



-ira 



{n\x\i)\ 2 {A* t (io ni )A t (io ni )). 



(5.13) 



Substituting the result of equation (5.11) into the above, we have 

o 

m 



$\ 2 = W \(n\x\i)\ 2 S(u; ni )t 



(5.14) 



From this, dividing the transition probability by time gives the transition rate 

,2 



Pi- 



t 



|c«| 2 = ^l(nkl^)r^Ki). 



(5.15) 



Therefore from equations (5.71 and (5.15) the transition rate becomes 



Tft ft 

p5(w ni )— [{i + l)5n ti+1 + i<5 n ,i-i] 
mS(uj r , 



2Huj 



n. 



i-l 



(5.16) 
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This gives rise to the total absorption power for the initial state \i) given by 



^ = ^[(. + Wl- ¥ [« 

that is, 

(5. IT) 

since S'(w) = £(— w) is an even function. Here we have dropped the subsrcipt of P since 
this absorption power turns out to be independent of the initial state. 

From this point we conclude the theoretical investigation which was carried out into 
the effect of a stationary time domain signal on a one-dimensional model of an atom 
trap. The derived relation is both concise and general, in terms of the power spectral 
density of the exciting force. The power spectral density term encodes full details of 
a broadband excitation containing a stationary signal as well as a possible additional 
stochastic noise. Now for the same spectral density in (5.17) when considering a simple 
harmonic oscillator, can be given as 

-it A 2 

S(u) = [S(w - 0) + 5{u + ft)} (5.18) 

where ft is the frequency induced by the electro-acoustic vibration generator, and to the 
internal trap frequency (i.e. the resultant frequency due to the combination of different 
magnetic fields from TOP and anti-helmholtz coils, A related to the acceleration of the 
shaker by A = —aft 2 . Hence, the absorption power can be expressed as 

V= 7 ^S(co-Q). (5.19) 

To obtain a definite yet simple estimate, we shall average P over a uniformly spread 
statistical distribution of the harmonic oscillators with density p(ui) = 1/0 for \u> — fl\ < 
SI /2 having a unity total probability. This yields the final expected value for the absorbed 
power by the atom trap as follows: 



/•Q+Q/2 

(V) = / V(oj)p(oj)duj 

JQ-n/2 

nma2Q3 (5.20) 



2 n 3 



As it can be noticed, the power produced by an atom is quite large on a timescale of 
a second compared to k B T for a temperature of micro Kelvin and to the energy of the 
ground state hf: 

k b T - 1CT 29 J 
hf - 1(T 31 J. 
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a 

Figure 3. Double-logarithimic plot of power (W) as function of a (mm). 

This indicates that the atom trap is extremely sensitive to an extremely small pertur- 
bation; Hence, even for an extremely small perturbation the power generated would be 
enough to produce evaporation within a few seconds. 

6. Discussion 

From the above description, it is clear that we could indeed find the BEC and the atom 
cloud reacting vigorously to the harmonic perturbation induced by the electro-acoustic 
generator. The brief analyse in the density of the Bose-Einstein condensate in section 4, 
mentioned the need for a more profound understanding of this new way of exposing the 
condensate. We found in section 5, that by deriving from a simple harmonic oscillator, 
the spectral density, and transition rate of the atoms in the trap can lead to determine 
the total absorption power of the initial state; and finally, using the same spectral density 
we realized that the power induced on an atom can be extremely large as shown in fig. (3) . 
This findings are good motivation for pursuing the experiment. 

On the other hand, there are numerous practical issues that will need to be resolved 
for the experiment to be successful. One of the problem is the influence of the trapping 
orbital electromagnetic field, which can potentially alter the dynamic of the condensate. 
For instance, we noticed that such configuration of the Quadrupole-Ioffe magnetic trap, 
when under the time-dependent potential of the electro-acoustic vibrator, will force the 
Majorana Hole to form a helix of the rotating TOP magnetic field. We need to know 
what combined effect it will have on the cloud or BEC. 

For most of the situation, s-wave length a s is not sensitive to the external potential, 
however, near Feshbach resonance the condensate can switch signs (i.e. repulsive to at- 
tractive), and change the dynamics of the condensate. The collapse of BEC can happen 
if the critical particle number N cr — N , that is, if during the oscillation a s is tuned by 
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the trapping frequency. Another problem may arrive due to mechanical vibrations of the 
equipments in the movable frame in the schematic. 

When it comes to flexibilities, one important factor is the way the increases on the 
energy level due to gravitational lamb shift will be measure; by using lasers off resonance, 
and reading the refractive index of BEC, which will differ by the presence of different 
energy levels, we will be able to tell to what extent the energy was increased. However, 
One of the constraints in this experiment is the number of atoms in the condensate nec- 
essary for the observation of gravitational Lamb shift; fO 6 atoms in the condensate of 
Rb87 may be slightly difficult to obtain and control for a long time. 

Finally, it is worth mention the difficulties that the experiment poses to verify the 
prediction of gravitational Lamb shift. However, difficult, the experiment is also viable, 
as long as the inertial effect will be dealt within the short lifetime of the BEC. Therefore, 
we feel the need to take the steps towards the development of the experiment, and study 
the effects to ensure that such discrepancies between gravity and quantum mechanics 
can be quantified and understood. And in addition to that, have experimental evidence 
that a clear observation of gravitational Lamb shift was achieved. 
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